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Abstract 

In the first part to this paper [l| it was shown how a simple Magnetohydrodynamic model could 
be used to determine the stability of a Tokamak plasma's edge to a Peeling (External Kink) mode. 
Stability was found to be determined by the value of A', a normalised measure of the discontinuity 
in the radial derivative of the radial perturbation to the magnetic field at the plasma-vacuum 
interface. Here we calculate A', but in a way that avoids the numerical divergences that can 
arise near a separatrice's X-point. This is accomplished by showing how the method of conformal 
transformations may be generalised to allow their application to systems with a non-zero boundary 
condition, and using the technique to obtain analytic expressions for both the vacuum energy and 
A'. A conformal transformation is used again to obtain an equilibrium vacuum field surrounding 
a plasma with a separatrix. This allows the subsequent evaluation of the vacuum energy and A'. 
For a plasma-vacuum boundary that approximates a separatrix, the growth rate 7 normalised by 
the Aflven frequency is then found to have ln(7/7^) = —^\n{q'/q). Consequences for Peeling 
mode stability are discussed. 
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I. INTRODUCTION 



Modern Tokamaks and future designs for power-plant scale Tokamaks, have plasmas 
with a separatrix at the plasma's edge. To determine whether these Tokamak plasmas are 
susceptible to an ideal Magnetohydrodynamic Peeling mode (External kink) instability, a 
simple model was generalised from a cylindrical to toroidal Tokamak geometry in the first 
part|l| to this paper. The conclusions from part (I) were: (a.) Peeling-mode stability is 
determined by the value of A', a poloidally averaged measure of the discontinuity in the 
radial derivative of the perturbation to the magnetic field at the plasma-vacuum interface 
(the separatrix), and (b.) that regardless of the sign of 6W the growth rate could still 
be vanishingly small. To determine the consequences of part (I) for the stability of the 
Peeling mode, this paper evaluates A'. It is essential that in this calculation any divergences 
due to the X-point are incorporated and not accidentally removed or accentuated by the 
discretisation of space that is usually required by numerical modelling, and ideally that 
A' is calculated exactly. That is the purpose of this paper. This paper and part (I) are 
summarised in Ref. j^. 

By using analytical methods to study Peeling mode stability in a plasma equilibrium with 
a separatrix boundary (and an X-point), we hope to avoid difficulties that are encountered 
with numerical studies, and to gain a better understanding of the physical factors that affect 
the plasma's stability. The techniques developed here may also find further applications to 
related problems using different models of the plasma. The results from this and future 
work are intended to provide understanding and tests, that will assist the development of 
codes for stability calculations that incorporate more advanced models of plasma physics 
and Tokamak geometry. It is also hoped that the methods developed in this paper will have 
applications outside of plasma physics. 

Outline: We review conformal transformations in Section [Tll and describe the Karman- 
Trefftz transformation!^ in Section lllli The Karman-Trefftz transformation provides an 
example of a transformation from a circular boundary to a separatrix boundary with an 
X-point. Section [IV] reviews how a complex potential may be defined and used to calcu- 
late how the vacuum magnetic field will transform under a mapping from a system with a 
circular boundary, to one with a separatrix. For a large aspect ratio system, the vacuum 
energy is calculated for both a circular cross-section and a separatrix cross-section in Section 



2 



\V\ obtaining the vacuum energy for a separatrix cross-section in terms of a sum of Fourier 
co-efficients. The Fourier co-efficients are determined by the plasma-vacuum boundary con- 
ditions, and Section IVTl discusses how these boundary conditions are modified by a conformal 
transformation. This is where we have departed from the conventional textbook applications 
of conformal transformations, that require the boundary conditions to be zero. Instead the 
transformed boundary conditions presented in Section |Vl] provide analytic expressions to 
determine the Fourier co-efficients in terms of the straight field line angle, that is not yet 
known. The straight field line angle is calculated in terms of an equilibrium vacuum field 
and the conformal transformation, in Section IVIII Section IVIIII calculates an equilibrium 
vacuum field for both the circular boundary and the separatrix boundary systems, subse- 
quently allowing analytic expressions for the safety factor q and the straight field line angle 
to be obtained at the plasma-vacuum boundary. At this point all the analytic expressions 
needed to calculate the vacuum energy have been obtained. Section [IX] investigates how 
other quantities (in addition to the field and the boundary conditions), change under a 
conformal transformation. These expressions are tested in Section |Xl by re-calculating the 
vacuum energy from a surface integral representation given in the first part to this paper. 
It is noted that this calculation infers the value for A' in terms of the vacuum energy 6Wv- 
Reassured with the results from Section |Xl Section |Xl] calculates A' directly, and confirms 
the same answer in terms of 6Wv found in the previous section. At this stage we have 
analytic expressions for SWy and A' in terms of a sum of Fourier coefficients, and analytic 
expressions for the Fourier coefficients. Section IXIII calculates 6Wv by evaluating the sum 
of Fourier coefficients in a number of ways, finding the same result as for an equivalent cal- 
culation in a circular cross-section system, and independent of the details of the separatrix 
geometry. Section IXIIII provides a discussion that compares and extends previous work, 
considers the mode structure, and summarises the paper. 



II. CONFORMAL TRANSFORMATIONS 

A conformal transformation w{z) (e.g. see Ref. j4|), is an analytic function between 
complex planes z —>■ w{z). It has the property that the angle (and direction of the angle), 
between curves in the plane from which they are mapped, is retained between the resulting 
curves in the plane onto which it is mapped. This angle-preserving property ensures that 
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the unit normal to a line in one plane will map to a vector normal to the mapped line, with 
a consequence that a boundary condition of = n^-Bz will be transformed to a boundary 
condition of = n^-Bw for the boundary in the w{z) plane. More generally, if the arc length 
around a boundary contour in the 2;-plane is parameterised by a, then it will be shown in 
Section EH that n^.V^V{w)\^ = |^| n,.W,V{z)\^. 

An important property of conformal maps is that a function that satisfies Laplace's equa- 
tion, will continue to do so after a conformal transformation. In other words, if V^V^(z) = 0, 
then provided w{z) is a conformal transformation, V^V"(2;(w)) = also. Riemann's mapping 
theorem indicates the existence of a conformal transformation from a circle to any closed 
region. So provided a suitable transformation may be found, and provided the boundary 
conditions map in a simple enough way (as they often do), then it is possible to find solutions 
in complicated geometries by solving a problem with a simple circular boundary. 



III. KARMAN-TREFFTZ TRANSFORMATION 



A mapping that may be used to take us from a circle to a shaped cross section with an 
X-point, is the Karman-Trefftz transformations!. The Karman-Trefftz transformation is a 
generalisation of the Joukowski transformation that is well known for its use in aerodynamics 
calculations for the lift from an airplane wing. It maps from a domain that surrounds a circle 
containing the point z = —I and whose edge passes through z = Z, to w(z), with 

' w + nl\ f z + 



w — nl J \z — I 

For simplicity we will restrict ourselves to domains of z that are symmetric about the real 
axis, so that w{z) is also symmetric about the real axis with a boundary that has 

z = -a + {a + l)e''^ (2) 

with a > 0, so that a will parameterise the boundary curve (in both the z and w{z) planes). 
Note that a is not the argument of either z or w{z) (if we are centred on z = a in the z-plane, 
a is the poloidal angle). If ri = 3/2, then the cusp-like point at z = I becomes an X-point 
(with a 7r/2 interior angle at the joining surfaces in the w{z) plane) [sl, and n = 2 produces 
the Joukowski transformation. By making l/a arbitrarily small we make the X-point region 



arbitrarily localised, a situation similar to that described by Websterjsj. This may be seen 
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by rearranging Eq. [T]to give 



and writing it as an asymptotic expansion in l/z, with 

So provided \l/ z\ is sufficiently small then w{z) = z. 

In summary, the Karman-Trefftz transformation provides an explicit representation of a 
transformation from a circular boundary (with z = —a + (a + 1)6^"), to a shaped boundary 
with an X-point. 

IV. THE COMPLEX POTENTIAL 

We will need to know how the vacuum magnetic field (the gradient of a potential) is 
transformed as we move from a circular cross-section to the X-point geometry. This is most 
easily accomplished by representing the magnetic field as a complex number whose real 
and imaginary parts are interpreted as its x and y components, and by defining a complex 
potential Q in terms of the magnetic potential V. 

The complex representation for the magnetic field is given in terms of the magnetic 
potential V, with 

B - — + (5) 

^ dx dy 

The complex potential Q is defined in terms of V and the conjugate function of V, such that 
Q is analytic and satisfies the Cauchy-Riemann equations. Specifically, 

n{z) =V{z) +iij{z) (6) 

with ip the conjugate function of V. Then the Cauchy-Riemann conditions are satisfied, 
with 

dV _ dip 

dx dy /y\ 

av ^ _di>_ ^ ' 

dy dx 

The Cauchy-Riemann conditions may be used to show that 



and hence the field in the transformed system w{z) may now easily be found from 



^ dQ dfl dz ^ dz , . 

^ dw dz dw ^ dw 



V. THE VACUUM ENERGY 

Working in terms of the complex magnetic field and the complex potential, we have 

SWv = / dw^dwy (10) 



where the integral extends from the boundary that is parameterised by a at w{z{a)), to 



infinity. To evaluate the integral we use Eq. [9] so that | p = | ^ f | ^ | ^ and we change 
coordinates back to the circular cross-section coordinates^ with dw^dwy = ^^^^^f^dxdy 
where ^''g^l'^^^ ~ because w{z) is an analytic function[4]. So when we change into the 
z coordinates (for the purpose of evaluating the integral Eq. [TD]), the factors of |^|^ and 
1^1^ cancel to give 

SWv = J \B,fdxdy (11) 

In a similar way it may be shown that Byj.dlyj = §^ Bz.dlz, refiecting the fact that the 
same total current is contained within 1^ and 

Hence the vacuum energy may be found in terms of the vacuum energy for a solution 
with a circular boundary, which is much easier to calculate. The actual values of 6Wv and 
Bz are determined by the plasma-vacuum boundary conditions, that will be modified by 
the transformation. The mapping of the boundary condition and the resulting boundary 
condition are obtained in Section IVIt but for the present we will obtain the general solution 
in terms of the Fourier coefficients that the boundary conditions will determine. 

The vacuum field has V A By = V.By = 0, so we may write By = W with V'^V = 0. 
We will start from a co-ordinate system with a circular cross section toroidal geometry, then 
subsequently obtain a 2D problem by taking the large aspect ratio limit. In this co-ordinate 
system 

VV = r \ \ (12) 

r dr dr da'^ B? dcj)'^ 

with r, a, the radial coordinate, poloidal and toroidal angle respectively. Writing 

p=oo 

V= e*^""*'"^^p(r) (13) 



and then projecting out the Fourier components, requires 

° = ;|:'-f^-'''^''W-"'(s)'^'W (14) 

which for any given (finite) n becomes 2-dimensional when r/R ^ 0, leaving 

nd'^Vr, dVn o . s 

° = 'V^ + "ar-^^- (15) 

that has solutions that tend to zero as r — > oo, olVp — ap (7^)'^', where the will be 
determined by the boundary conditions, and denotes the radial position of the plasma- 
vacuum surface. Note that it is only because the large aspect ratio limit makes the problem 
2-dimensional, that we are able to use conformal transformations in the calculation. 
The vacuum energy is 

5Wv^]^ J \B^\''dr (16) 
Into which we now substitute V, with 

J2 (^'^"'^% (— ) ^ (17) 
p=— 00 

and integrate with respect to and a, with r = at the surface, to get 

5Wv = 2Ti^R'^\p\\ap\^ (18) 

Here and in the remainder of this article R will be taken as a typical measure of the major 
radius that is approximately constant and independent of the poloidal angle. R is identical 
to the Rq of the first part to this paper, but because the rest of this paper considers a large 
aspect ratio limit, we will simply write R as opposed to it!o- 



VI. BOUNDARY CONDITIONS 

The boundary conditions in the ?iy'-planc determine n.B in terms of the plasma pertur- 
bation. We will write n.B in the w-plane as riyj.Byj. However, to obtain the co-efficients Op 
we need to know n.B in the 2;-plane (that we write as n^.B^). Therefore we need to know 

— * 

how n.B is transformed as we map between the z-plane where the boundary is circular with 
z{a) — —a + (a-|-Z)e"^, and the w-plane whose boundary is shaped and contains an X-point. 
This is calculated next. 
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Recall that the real and imaginary components are considered as orthogonal vector com- 
ponents. Then the tangent vector tw{a) of the surface traced by w{z{a)) is simply given 
by 

dw{z(a)) 

(19) 



t„(a) 



I 8^ I 
I da I 



However, using the fact that w{z) is an analytic function, so that = T^f^; then 



d'w(z{a)) 



dw dz{a) 



da 


dz{a) 


da 


1 dw 1 


dw 


1 — 


1 da \ 


dz{a) 


1 da \ 



dw 

dz 
I dw I 
I dz I 



(20) 



where the tangent of z(a) in the 2;-plane is again simply given by tz{a) = To 
obtain the unit normals we rotate the tangent vector by — 7r/2, by simply multiplying by 
—i. Hence 



We already know that 



w'{z) 
\w'{z)\ 

o obt 

w'{z) 
\w'{z)\ 



-it. 



w'{z) 
\w'{z)\' 



n. 



(21) 



Bz4^, so to obtain how riw-Bui transforms we need to simplify 



n. 



^ dw ) 



(22) 



where the dot product refers to the sum of the product of the real parts, plus the product of 

and ^ = r^^,/e*^»', to give 



the imaginary parts (examples may be found in Appendix IXVI) . We will use ^ = l/'^^ 

dw I I dw I 2 



dz 



f^/ Vf-\\ and write 



dz 

n,„.B. 



TBe" 



\w'{z)\ "-z-^z |«>'(z)|2 



cos(^^' + 9z) +i sm{9u,' + ^2)] . [cos(6l^' + 63) + i sin(6l^„/ + 63)] 
cos{6u,' + 9z) cos{6^i + 9b) + sin(6'^„/ + 9^) sin(6'^/ + 9b)\ 
l^cos{9z-9b) 

w' 

—Uz-Bz 

riz-Bz 



(23) 



I I 
I dz I 



This calculation is repeated by an alternative method in Section HXl 

Knowing how n.B transforms between z and the w{z) plane, we now return to the plasma- 
vacuum boundary conditions. As shown in part (I), the plasma- vacuum boundary condition 
is 



edge 



Vi).B 



V 



edge 



(24) 



(25) 

edge 



where "edge", refers to the equihbrium position of the surface. Because Bq.V^^ = Vip.Bi, 
we therefore require that 

edge eage 

For a single Fourier mode in straight field line co-ordinates = Cmi'ip)^^"^^ , with 6 = 
^J^udx', g = ^ § ^dx' , y = Jx the Jacobian of the orthogonal 4> co-ordinate 

system, with ifj the poloidal flux, x the poloidal angle, and the toroidal angle, and I{jp) is 
the flux function for which B = /V0 + V0 A Vip. After taking derivatives B.V^^ then gives 

Vij.B^ = {im - mq) ^^me*'"'"*"'^ (26) 



This may alternately be written as 

n^.i?.=zmAj^4^e^-^— ^ (27) 
RBp qR^ 



where A = 

Now we transform into the z coordinates, transforming both ny^.Byj = nz.Bz/\w'{z)\ and 
Bp = BpJ\w'{z)\, to get 

R \Bpz\ qR' ^ ' 

Using Eq. [13] along with 14 = (^k \^) ? and that n^.B^ = e^.VV^, then gives 

From which the Fourier coefficients are easily obtained by multiplying by and integrat- 



ing from a = — tt to vr, to give 



\\V\) R2nf''^ \Bpz\ qR^' ^^"^ 



In the following section we will see that 



which will allow us to integrate by parts once to get 



a„ = - 



I ^ime-ipa^^ (32) 



\p\J R27T 

To evaluate the coefficients Op, we will need an expression for 9{a), this is addressed in the 
following 2 Sections. 



VII. THE STRAIGHT FIELD-LINE ANGLE 



In the absence of equilibrium skin currents, the plasma's equilibrium field Bq equals the 



vacuum's equilibrium field Bq at the surface between the plasma and the vacuum, with 



B, 



edge 



edge 



. Therefore provided that we know the equilibrium vacuum field at the 



surface, then we also know the plasma's field at the surface. Consequently, if we know the 
vacuum field at the surface, then it is possible to calculate the straight field-line variable at 
the surface. Firstly we note that 



e 



udx 



qR^ 



Br, 



qR^ 



^ dl 

'b 



(33) 



p 



An element of arc length parallel to the tangent vector, dl^, has 



dl, 



dw 
da 



da 



dw dz 
dz da 



da 



dw 
dz 



dL, 



(34) 



Hence an element of arc length transforms such that \dl^\ = \^ \ \dlz\ = \^ \ rada, for 



a circular cross section of radius in the z-plane. Using this plus 
may write Eq. [33] as 



9{a) 



qR^ 



\dl^ 



\B. 



pw I 



qR^ 



dw 



dz 



da 



\B. 



\B,\/\w'iz)\, we 



(35) 



Hence if we know the equilibrium field, then we can obtain an analytical expression for the 
straight field-line coordinate as a function of a in the z-pla.ne. 



VIII. EQUILIBRIUM VACUUM FIELD 

The equilibrium vacuum field must (i) have a potential that satisfies Laplace's equation 
in the vacuum region, (ii) have n.Bo = at the plasma- vacuum boundary (including at the 
strongly shaped X-point containing equilibrium), and (iii) have the field -Bq = at the X- 
point. The first part is most easily satisfied - we can take a solution that satisfies Laplace's 
equation and n^-Bz = for a circular cross section, and after a conformal transformation to 
a shaped cross-section we will still have n^.Byj = and a potential that satisfies Laplace's 
equation. To obtain a field with Bp = at the X-point, we follow a procedure that is equiv- 
alent to that when applying the Kutta condition to obtain the fiow around an airplane wing 
(using a conformal transformation). Essentially, in the z-plane we combine a homogeneous 
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horizontal field and a circulating field, such that the field becomes zero at a single point on 
the circular boundary. This is physically equivalent to imposing an external horizontal field, 
and then driving a current through the plasma. Mathematically it corresponds to taking a 
complex potential in the 2;-plane of 

n = Bpo \i{z + a) - i ^ ^\ - 2i{a + I) \n{z + a)\ (36) 
I [z + a) J 

with a boundary at 

z = -a + {a + /)e^° (37) 

with a G [0,27r], BpQ a dimensional constant, and the radius of the circular boundary 
Ta = (a + /). The sign of BpQ determines the direction of the circulation of Bp^, clockwise 
{BpQ < 0) or anticlockwise {BpQ > 0), and the field Bp^ is obtained from Bp^ = with 



^^'-Tz- ' (7T^^^° ^^^^ 



which at the separatrix given by Eq. [371 gives \Bpz\ = 2BpQ{l — cos(q;)). To consider an 
outermost flux surface that is just inside the separatrix, we may instead consider 

z = -a+ {a + l- e)e^" (39) 

with e <^ /. Then for e ^ / ^ a we have 

\Bp,\ - 2Bpo (^1 - cos(a) + (40) 

so instead of Bpz = at the X-point (that is located at e = and a = 0), we have 
Bpz = /a?. Notice that we have retained the singular perturbation in e (singular in that 
although formally e^/a? <^ e/a, it is the term in e^/a? that qualitatively alters \Bpz\ by 
preventing it from being zero), further details are given in Appendix IXIVI The field in the 
transformed space is given by Bp^, = Bp^-^, although we shall not need this here. Plots of 
the equilibrium are given in Figure [TJ 

As we approach the separatrix the behaviour of 6{a) is dominated by the zeros in w'{z) 
and Bpz that occur near the X-point. Near the X-point it may be shown (in appendix IXIVj) . 
that for the case of n = 3/2 with field lines crossing perpendicularly to each other, that 
|w'(z)P is given by 
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FIG. 1: The figure shows contour plots of the imaginary part of the complex potential Eq. [36] for 
the magnetic field, with a = I = 1. Such plots give streamlines of the magnetic field J]. The plot 
on the left is equivalent to a combination of a vertical field and that produced by a current through 
z = i, with a field strength such that Bz = at the bottom of the plasma vacuum surface. The plot 
on the right is obtained from Q{z) by a conformal transformation, writing ^}{z{w)) and calculating 
plots of Q{z{w)) in the tn-plane. z{w) is obtained from Eq. [T]with n = 3/2 and / = 1, that gives 
z{w) = (^!|!3/2)2/3-(w-3/2)^/^ • plots the domain of z and w were rotated by substituting 

with iz and iw respectively, so that the X-point is seen at the bottom of the figure. 



We may use Eqs. HOl and |4T] to calculate q, with 

^ 271 27cR^J \Bpz\ " 27r R^B^o / ^1 - cos a + /2a^ ^ ' 

and y/c^ = 2V2 {lY {£). Similarly for e{a) we have from Eqs. [351 HOI and SH that 

eia) = 14^^ r , (43) 
q R^Bpo ^ v^l-cosa + eV2a2 



Because the integral is dominated by the divergence at the X-point where a = 0, q may be 
approximated by 

Ira I I da _ Civ^^inf^") (44) 



27r ^Jc^ J lo? I <?■ 2tx 

' ' 2 + 2a2 



with Ct = J/"^ ^ and similarly 



da Qv2 / e 
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Therefore as we approach the separatrix, with e — 0, g has a logarithmic divergence with 
q — ln(e) that is typical for a Tokamak plasma near the separatrix. Hence the qualitative 
features of 9{a) could have reasonably been postulated without showing that they also arise 
from a vacuum field whose potential satisfies Laplace's equation, has n^.Bpyj = on the 
separatrix, and with Bp^ ^ at the X-point; but it is reassuring to know that this is also 
the case. 

It is interesting to note that in this model for the equilibrium field, the angle at which the 
field lines meet at the X-point determines how strongly the divergence is there. For example, 
if instead of meeting at 7r/2 the lines make a cusp (tending to parallel as they meet), then q 
is finite (A cusp is obtained by taking n = 2 in the Karman-Trefftz transformation, Eq. [T]). 



IX. MORE TRANSFORMED QUANTITIES 



Now we start to return to our problem of calculating A', by firstly calculating 6Wv from 
its surface integral representation that is given in part (I), with 



TT 



52 



(46) 



where the integral is over the plasma surface. This requires us to know how n.V transforms. 
The calculation is done partly to reassure us that the transformed quantities are correct, 
but also because it is a simple step to subsequently obtain A'. 
First we calculate how V„, transforms. We have 



and 



V^/(^(w)) 



d d 



f{z{w)) 



(47) 



(48) 



where z = x + iy and w = Wx + iwy. Now we will use the chain rule to expand Vzf{z{w)), 
noting that because w{z) is an analytic function it satisfies the Cauchy-Riemann equations. 



dx 

dwx _ 
dy 



dWy 



dWy 



(49) 



dx 
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and in addition, = This gives, 



\ \ II owx ax owy ax aw^ ay awy ay 



dw^ / df I „• df 
dx \ dwx dwy 



dw. 



I awy_ 
dx 

■ dWy 



I 9f I _df. 



dx \ dwx dwy J dx y dwx dwy 

( dx 
dw 



(50) 



'I J \ dwx dwy 



dw y7 f 9^X7 f 

dx ^ wj — ^ -^J 



dx 

Hence we have — j^^^- 

Now we consider the transformation of ra.V. In calculating how more complicated ex- 
pressions transform, the author has found the following identities useful, whose derivations 
are given in Appendix IXVi 

a.bc = ab.c = db.c (51) 
ab.cd = {a.c){b.d) + {ia.c){b.id) (52) 
For example, to calculate n^.B^ we use Eq. [52] along with the results of Section |VT] to give 



n.B = Uw-Bo, 



w'{z) 



\w'{z) 

1 



\w'(z)\ 



dw " 



(53) 



- \w'(,z)\ 

as before. In the last step we used, 



dw dz 

dz ' dw 

dw ^ dz g 

dz ' dw 



that may easily be confirmed by writing ^ = a + ip , so that |^ = l/{a + iP), and multi 



(54) 



dz 



plying out. A similar calculation for n.V gives 
n.V = n,„.V„, = ^n^.Jv. 



1 

\w'{z)\ 



(55) 



X. RECALCULATING 6Wv 



Firstly we re-express Eq. |16] using our transformed quantities, then we show this gives us 
the same result, Eq. [18] from Section [V], before showing how the calculation easily generalises 
to give us A' in terms of SWy- 
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Using 



gives us 



V 



edge 



edge 



Vi).Bi 



edge 



edge 



This is substituted into Eq. |46l to give 

f m \ f I 



qn^ I \ m 



6Wy = 7^A<pj^dx{E]{l^.] Cme-'"''^''''' 



where A = Then using |. = we get 



\nq J yR^B"^ 
JL = "-V 



^ noip 



(56) 



(57) 



(58) 



5Wv = TTA{^y-f§dl (^) e-^-^+^«'^^n.V [RBpn.B, 



(59) 



where we used dl = J^Bpdx- Now we will transform this equation into coordinates in which 
the plasma has a circular cross-section, using: 





w'{z) \ dlz 


1 Bp^ 1 - 


|«,'(^)| 


Ti^.B^ 


- \w'(z)\ 




- \w'(z)\ 



(60) 



After using the chain rule to expand the term in n^.V^, we have 



(61) 



The Uz-V z operator acting on Uz-Bz will produce a term of order n larger than Uz-Bz- Thus 
usually we would neglect the second term. However here we need to be careful that there 
are no geometrically driven divergences, this is done in Appendix IXVIt where it is confirmed 
that the term is of order ^ smaller and may be neglected. Hence if we retain only the leading 
order term in n, then rearranging the expression slightly we have, 

5Wy = vrA f-V C /rfae--^(-)+^"^^^^^n,.V. (n..S,) (62) 
\nqj J B^m\w'{z)\ 

After comparison with Eq. [35] for 6{a), this may be written as 



5Wv = ttA 



nq 



2 t2 



\R-^B^ 



da- 



^—im0{a)+in4> 

imO'ia) 



-Uz.Vz {Uz-Bz) 



(63) 
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Using Eq. [HI 



^ipa—intf) 



a. 



(64) 



Giving 



'P ^2 



imO'ia] 



(65) 



Next we observe that 



da 



imO'ia) 



^ipa—im8{a) 



ip 



n 



(66) 



that we will justify below, and that appears to be the key result linking the high n and q 
calculations at arbitrary cross-section, to the circular cross section result. Integrating by 
parts we get 



da- 



da- 



^ipa—im8(a) ^ 



+ — (f> da 



0"(a) 



^ipa—imd{a) 



(67) 



im6'{a) j --^ ' ip J im{6'{a)Y 

To estimate the second term we notice that je^pa-^'^^C")! < i. Then taking 6{a) as given by 
Eq. HUthen we find 



X ^"(q) „wa-ime(a) 



< 



— f 

met J 



2 """2^ 



--da 



- -^V2(^) |ln( 

met \a/ \ \ 

„, (ilKlM 



27ra 



(68) 



where the integral is easily obtained by substituting ^ = Sinh(u). Hence for e/a ~ 1 the 
term is of order 1/n and may be neglected, and as e/a —> the term also tends to zero, and 
hence may be neglected. Thus in the high-n limit we get 

SWy = nA^-j e.(^^Ji^a, fdae (69) 



Taking the complex conjugate of Eq. [321 and rearranging, we get 



" P im^ 



(70) 



which upon substitution into Eq. [69l gives 



(71) 
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For the high m,n hmit considered here, we expect the |app coefficients to be largest for 
p m ^ nq, and hence in the high m, n hmit we expect 

5Wv = 2Ti^RY.,^^\p\\aA^ (72) 

where we have also taken P / {B?B'^) ^ 1. Hence we have re-obtained Eq. [T8]from the high-n 
expression given in part (I). This gives us confidence in the reliability of the calculations. 
In addition [|n2.V2(n2.-B^)|] may be estimated by approximating the plasma as a vacuum 
and solving Laplace's equation to approximate and obtain the perturbed field both inside 
and outside the plasma respectively, and correctly matching the fields at the plasma-vacuum 
boundary. Then we find that [|?T,2.V^(n2.i?;2)|] = 2nz.'Vz{nz-Bz). Hence the above calculation 
may be used to infer that the term — vrl^mpA^A' appearing in Eq. 70 of the first part to 
this paper, is equal to 26Wv, where 

A' = 

Hence to evaluate A', we need solely evaluate SWy- 



1 f P tttV^.^i 

^ i dlRB. i 
27r/ "R-'B^ vV-.^i 



(73) 



XI. DIRECTLY CALCULATING A' 



Here we show how to calculate A' directly, using the same assumptions as in Section |Xl 
For simplicity in all that follows we will take P/R^B^ ~ 1, and firstly use Eq. [2U that 



ViIj.B\ 



edge 



Vij.B^ 



to write 



edge 



A' 



dip 



Vtp.Bi 



— IdlRBp 

271 J " 



(74) 



Next we use Eqs. [601 to obtain 



A' = -3- i Tada — \— fln^.V in,,.B^, 



ZTT 7 Uz-Bz 



(75) 



Making the usual approximation that treats the perturbed field near the plasma's edge as 
behaving the same as in a vacuum, and also using Eq. [TTl gives 



Wz-^z [nz-Bz] 



dP 



e'^"'^^b|(b| + l) 



(76) 
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We also have 

n R — ^ I — — pi'pa-in4> \p\_ 

lI'z-J-'z dr\r=ra I^Pt^O^ "p r„ 

= Ep^o ^'""h § ime'{j3)e'^'^P^-'PPdl3 (77) 

where we used Eq. EJ] that imphes 9'{a) = Eq. M, and that ^'(a)e^™^(") = 

'^p^qG^'^"'-^ § 6'{(3)e^'^^^f^^~''P^dl3. Note that the above equation can be obtained more di- 
rectly from Eq. [28l 

Using the above results, after some cancellations we obtain 



Then using the result Eq. [66] of Section |Xl that 



we obtain the first term as 

_ 2 ^ B? ^ „ 1^1 {-iy) A 1 



,r$^A Z^P7^0"pl^'l J'"" ime'(a) " |5„|2 A2 Z^p^O "pl^*! |p| ^ 27r J ^ 

Using Eq. [32] for Op, and the same approximations as above, the second term gives 



(Later we will find that Ep^o bIkpP ~ m ^^"^2^ and hence that A' ~ —2m.) 



(78) 



^ da^-^e'P'- =i^§ e-.me+^pa _^ q (^^%^^ (79) 



(80) 



(81) 



^^^A Z^PT^O l/'l"? J" ime'ia) ~ ttUA ^Pt^O l/^ISip f ^ + ^ l^„g^ 

= Ep^o § e™^^'^)-*^^^/?^ / e'™^(")-*P"(ia 

Hence using all of the above, and Eq. [18] for ^VTy, 

A' - -AAE\p\M- + 4-2 f^J!l7ll (1 + I (82, 



XII. EVALUATING THE SUM 

We have found that the vacuum energy SWy and A' are both determined from 
E^-00 |p|kpp5 that we may in principle evaluate using our analytical expression for Op. 
We do that here. 
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Firstly note that if \p\ is replaced with p in Eq. [181 then we may easily resum the series, 
because 



A2 



27r 


i?2 














27r 


i?2 



-ime(l3) 



Jpl3J_ 



2tt 



(83) 



where in going from lines 3 to 4 we integrated by parts, and in going from lines 4 to 5 
we note that ^'(/3)e™^(^) = Ep^o^*^^^ § 9' {a)e'"''^^'^h-'P'^da. We might expect the values 
of the coefficients to be peaked for values of p ~ m ~ ng ^ 1, and so it is likely that 



liP|opP — EI^-ooPI'^pP- This has been confirmed by calculating the 



^00 

jp=—oo 



sums using a saddle point approximation. (The details of the calculation are too long to be 
included here) 

It is instructive to recalculate the result using a simple model for 6{a), that encapsulates 
the fact that as we approach the separatrix (with g ^ oo and the local field line pitch 
V becoming increasingly peaked near the X-point), the function d{a) becomes increasingly 
similar to a step function. In this simple model we take g ~ |, 5 <^ 1, and z/ ~ g ~ | when 
a G (—5, 5), this leads to a very simple model for 9{a)^ with 

a e (-TT, -5) 

{^)2ix aei-6,6) (84) 
2n a E {6, tt) 

So that 

1 

1 a G (5, tt) 

Because 6 {a) is piecewise linear, it is easy to evaluate ^ ciae*™''^"^-*^", that gives 

sin(p5) 



9{a) 



)2n 



a G (— TT, —5) 
a G (-5,5) 



(85) 



■ P 

-t—A—-m- 



\p\ R p{p6 



rmT) 



(86) 



Hence 



SWv 



2vr2^Ep^o \P\K 



2tt'^R^^^I ^^"^^Ep^O \p\(p5-mn)2 
R ^ I Jl IpI 



sin^ (pS) 



|(p5— m7r)2 " Jl \p\(p5+imr)'^ 



dp 



(87) 



2zl: 

R 



l^mpA^m^- for m > 1 
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the same as was obtained previously by the saddle point approximation. To obtain this 
result we used 

f sinV) ._r WW (88) 
~[ \p\{pS — rmiy Ji \p\{p6 — rmiy m \ rn? J m \ m J 

and 

y ^^'(p^) r ^^^l!MLrfp = o f o f ^1 (89) 

\p\ip5 — mrry Ji \p\(pS + mny \ m? J \ m J 

where we also used m nq j. Notice that the integrals do not diverge at p6 — nm = 0, 
because sm{pS) = sm{p6 — mvr) = for p6 — rmr = 0, this would not be the case if m were 
not an integer. Hence not only do we find agreement with the calculation using the saddle 
point approximation, but we again find that 



oo 



|p||ap|^ — i> pIo-pP as m ~ ng ^ oo (90) 



p=— oo p=—oo 

Therefore both methods suggest that provided m ^ 1 and n ^ 1, then 



oo 



E NISP- E P|a.P = A^^^ (91) 

p=— oo p=— oo 

In addition notice that for m ~ rag ^ 1 the result of neither approximation methods involve 
6 (or e), suggesting that the result may be generic and independent of the detailed form of 
e{a). 

Returning to the calculation of A', Sections 1X1 and IXll showed that at leading order 



A' = -2(^l^| (92) 



and that using the above results gives 

A' = -2m (93) 

In addition, the work described above suggests that the result is generic for perturbations 
with n ^ 1, regardless of whether the plasma cross-section is circular, or shaped with a 
separatrix boundary that contains an X-point. 
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XIII. DISCUSSION 



A. Scope & Purpose of the Calculation 

In the first part to this paper we started from the simplest model used to study Peel- 
ing modes, that considers Peeling modes in a cylindrical plasma at marginal stability, then 
generalised it to a toroidal plasma. According to the model, the energy principle's 6W is de- 
termined by the value of A', that is a normalised measure of the jump in the gradient of the 
normal component of the perturbed magnetic field. In this second part we have restricted 
ourselves to systems for which the vacuum magnetic field may be treated as being approx- 
imately two dimensional, as is the case for a sufficiently large aspect ratio Tokamak. This 
allows us to use a conformal transformation in our calculations, and at high toroidal mode 
number we have obtained analytic expressions for the vacuum energy and A', whenever the 
plasma is perturbed by a radial displacement consisting of a Fourier mode in straight field 
line co-ordinates. These expressions remain valid for a plasma cross-section that approx- 
imates a separatrix with an X-point, and appear to be generic, independent of the exact 
form for 6{a). Because it is possible to do this analytically, there is the possibility of making 
similar analytic progress with other linear plasma instabilities whose plasma equilibria have 
a separatrix with an X-point. Such calculations can provide physical understanding and 
useful tests during the development of codes to study the stability of more general geometry 
Tokamak plasmas, either giving confidence in a code or indicating its limitations. 

B. (In) Stability of the Ideal MHD Peeling Mode? 

According to the model developed in part (I) and our calculation here of A' = —2m, we 
can now examine the model's predictions. According to the model developed in part (I), for 
the trial function used by Laval et al , stability is determined by the sign of 



2 



sw 



A AA' + J 



(94) 



R 



with 



A 



m — nq 



(95) 



nq 




(96) 
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A' 



1 /-...^ P M^^-Bi 



dlRB„ 



2nT—'R^B^ V^.fii ^^^^ 
and as mentioned in Section |Vl R = Rq is constant for the large aspect ratio limit considered 
here. 6W is minimised with respect to A (or equivalently, a particular choice of toroidal 
mode number), finding A = — J/(2A'). Then using our result of A' = —2m we get 

When checking the dimensions of 5W it is essential to remember that = iyi'-Om ~ rRBp, 
then because 5^ ~ p is an energy per unit volume, SW ~ r^RB^ and hence has units of 
energy. For the process of minimisation n was treated as a continuous variable (that for 
m ~ ng ^ 1 is a reasonable approximation). Now we consider two cases in turn, firstly 
V0. J = 0, for which 

I ri -I'Bl 1 /• -//' 



p 



^ f d^B,^ - 1 (99) 



2tx J RBp fi2 2n J ' RB"^ 

So that although 5W < for all m, because m ~ ng — oo then 5W ~ ^ ^ 0. For the case 
of V(f).J 7^ however, 

f 1 r I j.B 1 r di I ^j.B I r.J-B\ 1 r di I 

'=2^f''RB;^=2^f Bp-R^""^ " \^ 5^ /2^fB;R-^=' \^ 5^ 

(100) 

where () denotes the poloidal average, and q = ^ § ^-§- is the safety factor p. The 
divergence in at the X-point will mean that the poloidal location of the X-point will 
affect the value of J (as a function of g), this is also the case for Mercier stability and is 
discussed by Webster [5]. Then using Eqs. [981 1100[ and m ~ ng we have 

5W ~-(|)^Mi(i?4|)2<o (101) 

Therefore if 5W < is taken to indicate instability, the result would indicate that the 
Peeling mode remains unstable near a separatrix. However as observed in the first part of 
this paper the growth rate 7 has 7^ = —5W/ J (irpol^Pj with J c?rpo|^P diverging at a rate 
proportional to q'{ip). This gives ln(7) = — iln(g'/g) for the limit of a separatrix with g 
and g' tending to infinity, so that although SW is non-zero and negative, the mode will be 
marginally stable. This is similar to the calculation of the Mercier coefficient by Webster jsl, 
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where Dm is found from the ratio of two diverging quantities, with Dm ~ g/g' — as the 
separatrix is approached. Note that these resuhs for A' and the growth rate appear to be 
generic and independent of the detailed forms of 6{q.) or the radial structure of the mode. 

C. Is the Mode Physically Acceptable? 

Because we require a poloidal mode number m ~ ng, then near the separatrix where q 
oo we also require m —>■ oo. This raises the question: Is the mode physically acceptable? To 
answer this we reconsider the trial function, that has ^ ~ gime{a)^ with 9(a) = ^ i^{ce)da. 
When m nq the trial function becomes ^ ~ e™-^ from which we can see that because 

V is of order one and well behaved everywhere except near the x-point (where it diverges to 
infinity), then so also is the mode's structure. As discussed in part (1), the divergence in 
m ~ ng — > cxD is only manifested in close proximity to the x-point where the divergence in 

V causes the mode to oscillate increasingly rapidly as the x-point is approached. Elsewhere 

V is typically of order 1, and the mode structure is like that for a finite mode number, 
oscillating at a modest rate of order nv -C ng, and only weakly affected by the proximity of 
the flux surface to the separatrix. Hence the mode has a simple structure everywhere except 
for a region close to the x-point where it oscillates so rapidly that MHD would no longer 
be applicable. The resulting mode structure is consistent with observations of ELMs 71] , 
that show filamentary structures that follow the magnetic field lines, and whose poloidal 
structure near the X-point is difficult to determine. 

D. Previous Analytical Work 

As mentioned at the outset, Laval et al 6| considered a trial function consisting of a single 
Fourier mode in a straight field line co-ordinate, that is resonant at a rational surface in 
the vacuum just outside the plasma's surface. For that trial function, they found that for 
a positive non-zero current at the plasma edge, 5W < 0, and suggested therefore that the 
Peeling mode would be unstable for a non-zero positive current at the plasma's edge. On 
the basis of the sign of 5W our study also finds this, but our study also suggests that the 
growth rate will asymptote to zero as the outermost fiux surface approximates a separatrix, 
so that the mode will be marginally stable. 
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Lortz'8] also considered Peeling mode stability, in toroidal plasmas with shaped cross- 
sections, using a systematic calculation with a trial function whose resonant surface is inside 
the plasma. An advantage of the calculation by Lortz[8], is that the radial structure of the 
mode is considered. An unfortunate complication for this discussion is that in the ordering 
scheme of Lortz[8|, the vacuum energy can be neglected. This was not the case for our 
calculation in Part or of Laval et al|^. Nonetheless, we will consider the predictions 

of this calculation in the limit where the outermost flux surface approximates a separatrix. 

Connor et alfsl review the calculation of Lortzjsl, and use it to consider trial functions 
with resonant surfaces both inside and outside the plasma. They find that stability of the 
Peeling mode requires 



1 - ADm > 2 



S 

I 

'P 



1 



where the Mercier co-efficient D 



M 



-Q/P, and P, Q, S are defined as. 



(102) 



P = 2niqr [§ m^dX 



X 



2n J R-^Bl 



(103) 



S = P + q' § 



§ 



R^Bl '^xdx 

Substituting Dm = —Q I P into Eq. 11021 gives the stability requirement 

2 



Q 



p 



s 



p 



s 



p 



(104) 



Substituting for P, Q, and 5", allows Eq. 11041 to be simplified to 

J- 



P 



^1 



iy.wz ^ Jxdx-P'iz ^ i^dx 



2ti 



> (105) 



This may be simplified further by noting that because u = IJ^/R^, and in a large aspect 
ratio ordering where R is taken as approximately constant. 



d 1 



d 1 



vR? 



dx 



~2^ J d^j ^ 



R%^ ^ vdx 



and 



The Grad-Shafranov equation in x, co-ordinates, has 



dv 

dip 



V 



d_ 

dip 



[p + B') 



R^B^ d 
I 'di) 
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R^ 



(106) 



(107) 
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Therefore if R is taken as approximately constant (as would be the case either in a large 
aspect ratio limit or if we are sufficiently close to the separatrix that the integral is dominated 
by the divergence at the X-point and R may be approximated by its value R = Rx there), 
then using Eq. 11081 Eq. 11051 simplifies to a condition for stability of 

This may be simplified further still, to 



(109) 



27r 



B^R^ 



2 V0.J - 



dil) 



dx 



(110) 



Because the integrals are dominated by the divergence oi u/B^ at the X-point and < 
at the X-point, then based on the formulation of Lortz 8|, |9|], then provided V0.J > 
the negative expression clearly indicates instability to the Peeling mode. However, if we 
allow V0. J to be negative, then the formulation of Lortz ^, 9] also suggests that stability is 
possible provided 



< 



2n 









952] 
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V0.J 


p 






dip 



dx 



;iii) 



Therefore in principle there is a range of negative current values at the plasma edge for 
which the Peeling mode is stable. The appendix calculates dB^/dip near the X-point for 



a standard and a "snowfiake" [13] divertor. Interestingly, whereas a conventional X-point 
has a range of negative current values for which the Peeling mode is stable, in the limit 
of an exact "snowfiake" X-point (with fiux surfaces meeting at an angle of vr/3), the range 
of values of negative current for which the Peeling mode is stable, tends to zero. Whether 
this observation will have consequences for the plasma behaviour in a "snowfiake" X-point 
geometry remains to be seen, but it is a qualitative difference between a conventional X-point 
and that produced with a "snowfiake" divertor. 

As mentioned previously, the calculation also considered the radial structure of the mode, 
with ^ ~ x'*'='=, X a radial co-ordinate, and 

1 



A 



± 



± 



As we approach the separatrix, Webster [5|] shows that D 



1 + ^ 

4 P 



(112) 



0, giving 



Ah 
A_ 



— p ^ u 



(113) 
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and mode structures of ~ ^ and ~ x^l^ . For the perturbations to satisfy the 
boundary condition of a mode amphtude that tends to zero in the plasma, this requires 
us to use ^_ for resonances outside the plasma (the "external" Peeling mode), and ^+ for 
resonances inside the plasma (the "internal" Peeling mode). It should be noted that a 
potential problem with the analysis of Lortz et al[8j] when applied to Peeling modes, that 
the mode is taken to be sufficiently localised that the equilibrium quantities (that include q 
and g'), are approximately constant. This is almost certainly not the case near a separatrix. 

E. Summary 

We have started from a simple model for the Peeling mode, at marginal stability in cylin- 
drical geometry, and in Part (I) of this paper generalised it to toroidal Tokamak geometry. 
A conclusion of Part (I) is that Peeling mode stability is determined by the value of A', 
a normalised measure of the discontinuity in the gradient of the normal component of the 
perturbed magnetic field at the plasma-vacuum boundary. Therefore this paper evaluated 
A' in a large aspect ratio Tokamak geometry with a separatrix and X-point, but in such 
a way that the effect of the X-point is captured exactly, without encountering the usual 
discretisation errors present in most numerical methods. This was possible by generalis- 
ing the method of conformal transformations beyond textbook presentations, that require 
a boundary condition of either the function or its normal derivative to be zero. Here we 
observe that even if the field's normal derivative is non-zero at the boundary, it is still possi- 
ble to use the conformal transformation method. In this case instead of obtaining an exact 
analytic solution (as would be the case if its normal derivative were zero on the boundary), 
the 2-dimensional problem is reduced to a 1-dimensional problem that may subsequently 
be solved exactly or approximated. The approach avoids the errors that may arise due to 
the discretisation of space near an X-point, that are necessarily present in most numerical 
methods. This paper also calculated analytical expressions for physically realistic examples 
of the equilibrium vacuum magnetic field, and the straight field line angle at the plasma- 
vacuum boundary. These and other results are likely to find opportunities for application 
elsewhere. 

It is found that a radial plasma perturbation consisting of a single Fourier mode in 
straight field line co-ordinates with a high toroidal mode number n, in a plasma equilibrium 
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with a separatrix and an x-point, will produce the same change in the vacuum energy as 
the equivalent perturbation in a cylindrical equilibrium, with 5Wv = 27r^J^^A^m. It also 
results in the same value for A', with A' = —2m, where m is the poloidal mode number. 
Despite our trial function requiring m ~ rig — oo, we observe that the trial function has 
^ ~ ^^^!^vdx i\^qx is physically well behaved for all but a highly localised region near the 
X-point where MHD will fail to apply. Therefore we believe the trial function is physically 
acceptable even for a separatrix boundary. 

Previous work by Lortz sl and Connor et al [o}] was considered for an outer flux surface 
that tends to a separatrix with an X-point. Like Laval et al 6(] their work predicts the Peeling 
mode to be unstable if there is a positive current at the plasma's edge, and it also finds a 
well behaved radial structure for the mode. Interestingly, for a conventional X-point there 
is predicted to be a range of small but negative edge-current for which the Peeling mode 
is stable, but in the limit of an exact snowfiake divertor this range shrinks to zero size - a 
qualitative difference between a conventional and a snowfiake divertor. A limitation of the 
Lortz calculation is that it approximates the equilibrium quantities as constant on the length 
scale of the plasma instability, this is not necessarily the case for q or q' near a separatrix, 
and therefore the results when applied to a separatrix case should be treated with caution. 
Likewise, as noted in Part (I), there are potential limitations to the high-n ordering form of 
5W used here, and this should be investigated in future work. 

Thus we have developed a simple model for the Peeling mode, and found that despite 
5W < 0, the growth rate 7 tends to zero as the outermost flux surface tends to a separatrix 
with an X-point. As the outermost flux surface approaches a separatrix, the growth rate 
falls with ln(7/7yi) = — iln(g'/g); this has subsequently been confirmed with ELITE (S. 
Saarelma, private communication), leading us to believe that the effect of a separatrix on 
the high toroidal mode number ideal MHD model is now understood. The ideal MHD 
prediction of marginal stability at the separatrix means that other non-ideal terms such 
as resistivity, non-linear terms, or terms neglected in the high-n analysis, will play a role 
in determining the eventual stability. In general it is hoped that the methods and results 
contained in this paper will provide new tools for studying plasmas in separatrix geometries, 
and have potential applications in future studies of plasma stability and more generally 
outside of plasma physics. 
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XIV. w'{z) NEAR THE X-POINT 

To obtain w'{z) we differentiate both sides of Eq. [T]with respect to z, and rearrange the 
resulting expression to get 

^n— 1 
{Z + if^ 

and hence 



w'{z) = {w{z)+nlY^^^l-, (114) 



\z-l 



2(n-l) 



w'{z)\' = \w{z) + nil"" (115) 



\z + 0''-''^ 



We have deliberately obtained an implicit expression for |w'(2;)|^, with w'{z) given in 
terms of w{z). This is because one cannot simply expand w'{z) in powers of e, because 
the expansion will give the incorrect answer as e ^ compared with the exact result for e 
small but non-zero (i.e. e 7^ is a "singular perturbation"). Instead by obtaining |w'(z)p 
implicitly in the form given by Eq. I115[ we need solely be careful with the term \z — 
because \w + nl\^ and \z + are well behaved when expanded in e, as e — 0. Near the 

X-point, 

\w(z)+nl\^ = (2n/)4 + 0(e) 
|^ + /|2("+i) = (2/)2{"+i) + 0(e) 
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We need to be more careful with |z — that with 



-a + (a + / — e)e*° gives 



\z-l\ 



((a + / — e) cos(q;) — a + Z)^ + (a + Z + e)^ sin^(a) 
2(a + l){a + / - e) 



cos(a) + 



2(a+l)(a+l-e) 



;ii7) 



2a2 



1 — cosfa) + 



2a2 



where we have retained the term that provides the singular perturbation that prevents (z—l)'^ 
becoming zero for e 7^ 0, but neglected all the lower order terms that modify the answer by 
of order e/a and I /a. In principle some plasma cross-sections might require the retention of 
terms of order I /a, but here we neglect them so as to keep algebraic details to a minimum. 
Therefore at leading order we have 



\w'(z)f ~ 



(n-1) 



(118) 



" (2/)2{«+i) 

that for the case we are most interested in here with n = 3/2 (corresponding to an X-point 
with a 7r/2 interior angle), we have 



\w'iz)t 



i3iy 



(2a2)'/' ( l-cos(a) + 



(2/)5 
which is Eq. SH 



cos(a) + 



2a2 
(119) 



XV. SOME IDENTITIES INVOLVING COMPLEX NUMBERS 

In the following we write a = + iay, b = + iby, c = + iCy, and d = + idy, and 
remind the reader that the dot product refers to the sum of, the product of the real parts 
plus the product of the imaginary parts. For example a.b = aj)^ + ayby. Multiplication is 
as usual for complex numbers, for example ab = axbx — CLyby + i{axby + ayb^)- Then we find 

ab.cd = [{a^bx - ayby) + i^a^by + ayb^)] ■ [{c^dx - Cydy) + i^c^dy + Cydx)] 

{^^x^x ^y^y^i^^x^x ^y^y^ ~l~ is^x^y ~l~ ^y^x^ij^x^y ~l~ ^y^x^ 

^x \px(,^xdx ^y^y^ ~^ by(^Cxdy H~ Cyrf-j;)] -|- Gy [bxi^Cxdy H~ Cydx) by(cxdx ^ydy^\ 

= [dx{bxCx + byCy) + dy{byCx - Cybx)] + ay [dy{bxCx + byCy) + dx{bxCy - byCx)] 

{axdx ~\~ aydy) (b^Cx ~\~ byCy) -\- (^axdy aydx)(byCx Cybx) 
= {a.d){b.c) + {ia.d){ic.b) 

(120) 
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and similarly for a.bc, 



a.bc = a. {{b^c^ - byCy) + i{b^Cy + byC^)) 

^x^x^x (^x^yCy H" ^y^x^y H~ OjybyCx 
C^xis^x^x H" ^y^y^ ~^ C^yi^O^ybx ^x^y^ 

= (c^ + iCy). {{a^rbx + dyby) + i{aybx - a^^by)) 
= (c^. + iCy). {a^{bx - iby) + ay{by + ihj)) 
= {cr, + iCy). {a^ibnc " iby) " iay{-br, + iby)) 
= c. {{a^ + iay)% - iby)) 
= cab = ab.c 



(121) 



XVI. 2ND TERM IS ORDER ^ SMALLER THAN SWy 

Here it is shown that the second term in Eq. [611 here written as 5Wg, is of order 1/m 
smaller than 6Wv- The term we are interested in is 

\Bpz\ 



Using Eq. Hi] and Eq. HO we get 



n 



(nz.B,)nz.V;, 



\w'{z)\' 



(122) 



n 9 



2V2 (1)^ {^1 - cos(a) + eV2a2 



A/l-cos{o)+e2/2a2 



Ira 



(123) 



where 
gives 



^ -| 2 

2-\/2 (I) and is a constant. Substituting the above Eq. 11231 into Eq. I122[ 



\nq J \ a"^ J n 



Integrating by parts then gives 



52 



\nq J \CL J nm J 



riz-Bz — 
da 



driz.Bz 



da 



Using Uz-Bz = ^ and Eq. [T7|we get 



_d_ 
da 



(nz-Bz] 



E 



-ip\p\a. 



(124) 



(125) 



(126) 
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Substituting this into 11251 we get 

SWo = >r A f ^) C ( ^) c.^^ E / e-'-H-"* (127) 



p/O 

Where the poloidal dependence in ^ has been neglected due to the large aspect ratio 
Tokamak ordering. Using Eq. [32] we get 



^+^v^da = -^a*JPi-^27^ (128) 



which may be inserted into 11271 and after some cancellations gives 

The sum X^p^o 1*^?!^ '^^y evaluated arbitrarily accurately, as is indicated next. We 
use the alternative expression for ap given by Eq. [301 of 

A U 1 



" \p\ R 2n 



^ ime\a)e'"'^-'P''da (130) 



which gives 



^p l/^l rpl " R2 Z^p^O 

= A^^fJ ^ e'{(3)e-'^'^^^ Ep^o e+'^^2^ / 6' {a)e'"''^"h-'P'' (131) 

= A2^f^^(e'(/3))^d/3 

where in the penultimate line, we resum the Fourier series by noting that 9'{P)e^"^^^^^ = 
X]p^o § e~*P°6''(a)e™^*^°^ We approximate § (6''(/3))^ using the analytic expression 

for 9{P) obtained from Eqs. [35], [iQl and [HI obtaining 

^ ^^^^ ^ " c. / 1 - cos(/3) + eV2a2 " c. i^^E^^o e ^'''^ 

Hence taking = a + / ~ a we get 

= 2vr^%^A^ = (^) 5iyv 

where the last line uses the expression for 6Wv ~ A^m obtained from the main text. Hence, 
we may neglect this term compared with SWy- 
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(133) 



XVII. MODIFICATIONS TO PEELING MODE STABILITY FOR A 
"SNOWFLAKE" DIVERTOR 



The integrals that determine the stabihty of Peehng modes, are dominated by the diver- 
gence in v / B"^ that occurs near the X-point in the separatrix. This allows the integrals to 
be estimated by expanding the poloidal flux functions in the vicinity of the X-point, as has 
been done in Ref. [121 for both a standard X-point and a "snowflake" divertor's X-point. 
Near a standard X-point, Eq. 15 of Ref. [12| gives 




(134) 



with B havin g th e poloidal field's dimensions and L a length scale, and the poloidal magnetic 
field given hy B^ = and B^ = leading to 



(136) 



Note that near the X-point, Vip — > i?xV$, with Rx the major radius at the X-point. 
Therefore near the X-point = and we calculate dBp/dtp from 



dBl 

d<S> 



(136) 



and use Eq. 11341 to substitute for = 2 (^^j + S-z^, giving dB^/d^ at constant $ 



as 



dB^ / B 



+ (^^L/B 



(137) 



and similarly Bp = \^B /L j 2 + ^L/Bj. The minimum value of 2; is at a: = 0, giving 
2 



z^ain = -^^L/B, and 



dB^ 



9$ 



B 



(138) 



For a "snowfiake" divertor, Eq. (2) of Ref. [I2j] gives 



AI 



X z — 



(139) 
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with I the plasma current, A has dimensions of the inverse cube of the length scale over 
which the poloidal magnetic field varies near an X-point, and c is the speed of light. In a 
similar way to the calculation for the standard X-point, this leads to 



9$ \ C J ) 



with z^-^ = -3^c/AI, giving 



9$ 



fAI\ 4 /c$\ 



that tends to zero as we approach a separatrix. Therefore, whereas a conventional X-point 
leads to a Peeling mode stability boundary for which the Peeling mode can be stable for a 
range of small but non-zero negative current at the plasma's edge, this window of stability 
tends to zero size for a snowflake divertor. 

Finally, for a "snowflake plus" divertor [l2], a similar calculation using the Eqs. of Ref. 



12] finds dBp/d^ ~ y "1^ with I the current in the divertor coils and I^o the divertor coil 
current required for an exact snowflake divertor. In that case the range of values of negative 
edge-current for which Peeling modes are stable, tends to zero as / — > Ido- 
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